Nonequilibrium stabilization of charge states in double quantum dots 
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We analyze the decoherence of charge states in double quantum dots due to cotunneling. The 
system is treated using the Bloch-Redfield generalized master equation for the Schrieffer- Wolff trans- 
formed Hamiltonian. We show that the decoherence, characterized through a relaxation r r and a 
dephasing time t^, can be controlled through the external voltage and that the optimum point, 
where these times are maximum, is not necessarily in equilibrium. We outline the mechanism of 
this nonequilibrium-induced enhancement of lifetime and coherence. We discuss the relevance of 
our results for recent charge qubit experiments. 
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The loss of quantum coherence is a central paradigm 
of modern physics. It not only governs the transition be- 
tween the quantum-mechanical and the classical world, 
but has recently also gained practical importance in the 
context of engineering quantum computing devices. De- 
coherence naturally occurs in small quantum systems 
coupled to macroscopic heat baths. A huge class of such 
baths generates Gaussian noise and can hence be mapped 
on an ensemble of harmonic oscillators as in the spin- 
boson modeli. This can even apply, if the fundamental 
degrees of freedom of the bath are fermions, as it is, e.g., 
the case if the bath is a linear electrical circui1»2*2i, which 
is producing Gaussian Johnson-Nyquist noise. 

In this Rapid Communication, we study a generically 
different system: a double quantum dot coupled to elec- 
tronic leads. Such systems are studied as realizations 
of quantum bits^&*&. The position (either left or right 
dot) of an additional spin-polarized electron is used as 
the computational basis of a charge qubit as realized in 
Refpi. For another proposal of a charge qubit in semi- 
conductors, see Refi 

Our system simultaneously couples to two distinct 
reservoirs of real fermions. Other than oscillator bath 
models, this allows for the application of a voltage be- 
tween these reservoirs as a new parameter for control- 
ling decoherence. The voltage creates nonequilibrium be- 
tween the baths, which to the best of our knowledge has 
not been studied yet in the literature on open quantum 
systems. 

We study the dynamics of the reduced density matrix 
and identify the usual two modes of decoherence, dephas- 
ing and relaxation: Dephasing is the loss of phase infor- 
mation, manifest as the decay of coherent oscillations. 
This corresponds to the time evolution of the off-diagonal 
elements of the reduced density matrix in the energy ba- 
sis. Relaxation is the process during which a quantum 
system exchanges energy with the environment and ends 
up in a stationary state. This is described through the 
evolution of the diagonal density matrix elements. We 
are going to show that, surprisingly, the charge states 
can be stabilized by external nonequilibrium, i.e., the re- 
laxation time is longest at a well-defined finite voltage. 
We will show, that this working point is also very favor- 



able in terms of dephasing but competes with another 
local maximum at zero voltage. Our theory should also 
have applications in other systems. 

We consider a double quantum dot system with an ap- 
preciable tunnel coupling between the dots allowing for 
coherent molecular states in these systems^. The compu- 
tational basis is formed by the position states of an ad- 
ditional spin-polarized electron^. A superposition can 
be created by variation of the interdot coupling. In order 
to stabilize the charge, the coupling of the dots to the 
two leads is driven to weak values and the dot is tuned 
to the Coulomb blockade regime^ where sequential tun- 
neling is suppressed through the addition energy. Even 
then, the system couples to the environment through 
cotunnelingiS, the correlated exchange of two electrons 
with the external leads which ends up in a state with the 
same total charge as the initial one. 

The relevant Hilbert space is spanned by four states 
written \i,j) denoting i,j additional electrons on the left 
and right dot, respectively. |1,0) and |0, 1) define the 
computational basis as they are energetically accessible, 
the closest virtual intermediate states for cotunneling are 
|wo) = 1 , 0) and \v2) = 1 1 , 1 ) - This model applies if all 
relevant energy scales of the system (e as and 7, see below) 
are much smaller than the charging energies to the next 
virtual levels (e V2 and e Vo , also below) , which in turn have 
to be smaller than the orbital excitation of the individual 
dots. This can be realized in small dots. 



The total Hamiltonian of this system can be writ- 



ten as H = Hq + Hi where Hq 



H ICS de- 



scribes the energy spectrum of the isolated double-dot 



through H s 
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a" + a^'a") and the two electronic leads H lcs = 
Ecer&r^r + E^r^r- Thc sum over the dot 

k k k *-<kl ki kl kl 

states only runs over the restricted Hilbert space de- 
scribed above, the a L ^ R act on the lowest additional 
electron state on either dot. The double-dot is char- 
acterized by the asymmetry energy e as = ei — e r be- 
tween the individual dots and the interdot tunnel cou- 
pling 7. The virtual states |t>2) and \vq) are separated 
from the system by energy differences e V2 (upper virtual 
level) and e Vo (lower virtual level). The tunneling part 
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describes the coupling of each dot to its lead and will 
be treated as a perturbation. t c represents the tunnel 
matrix elements between the dots and the leads. It can 
be absorbed in a tunneling rate hT = 2TTt 2 N(ep). This 
has to be chosen small such that the Kondo tempera- 
ture is low Tk <C T and perturbation theory holds, e.g., 
T = 10 9 Hz. Figure H] shows a sketch of the setup. 
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FIG. 1: (Color online) Sketch of the considered artificial 
molecule in the Coulomb Blockade regime, where 28 is the 
level splitting and V = [ir — \il the bias voltage, that is ap- 
plied between the two leads (hatched). The virtual states \v2) 
and |«o) are outside the plotted energy range. 

From now on, we use the basis of molecular states 
obtained by diagonalizing H sys with splitting 25 — 
2\J e 2 s + 7 2 . In order to treat cotunneling by leading- 
order perturbation theory, we rewrite H\ using a 
Schrieffer- Wolff transformation^. This removes the tran- 
sitions to the virtual states and generates an effec- 
tive Hamiltonian containing indirect transition terms be- 
tween the molecular states. A more detailed description 
of our calculation is given in Refs.H and El The final 
Hamiltonian is of the form H = Hq + H[ where 
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where Y and Y' denote right or left lead, the as describe 

Y Y' 

molecular states and the H-'- are given through 2nd 

k,k',c,d ° ° 

order perturbation theory, i.e., they are of 0(T 2 ). Note, 
that H[ conserves the particle number because it acts 
upon the double-dot by injecting and extracting an elec- 
tron in a single step. The terms with Y ^ Y' transfer 
charge between different reservoirs. Note that Eq. is 
a simple and generic Hamiltonian connecting a quantum 
system to two distinct particle reservoirs and is poten- 
tially relevant for systems other than quantum dots as 
well. 

We study the open system dynamics in the case of a 
time-independent Hamiltonian with a fully general ini- 
tial reduced density matrix. We use the well-established 
and controlled Bloch-RedfielcU^, which has been demon- 
strated to work down to low temperature for certain 
models 1 ™. It involves a Born approximation in H[, i.e., it 
captures all cotunneling processes in lowest nonvanishing 
order. The Redfield equations 1 - for the elements of the 



reduced density matrix p in the eigenstate basis of H sys 
(i.e., the molecular basis) read 

Pnm{t) = —iWnmPnm(t) — ^ RnmklRkl{t) , (2) 

k,l 

where u> nm — (E n — E m )/% and the Redfield tensor ele- 
ments R nm ki are composed of golden rule rates describing 
different cotunneling processes, which are essentially in- 
dependent due to the low symmetry of the system. Each 
contribution has a typical cotunneling structur o 5 ; 12 . An 
overview of the most important processes is given below. 
n, m, k and I can be either + (excited molecular state) or 
— (molecular ground state) with according energies E±. 
This type of perturbative analysis is only valid above the 
Kondo temperature Tr^, which can be easily driven to 
low values by pinching off the tunneling barrier to the 
leads. 

From the formal solution of Eq. (|2J) we can identify the 
relaxation and dephasing rates as 
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r r — Rc(i?+_|_+^ + R 

T = Re(i?+_ + _) = Re{R. 
The transition frequencies w nm are weakly shifted. 
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FIG. 2: Examples for relevant processes in the system: (a) 
a relaxation process that carries no current, (b) a relaxation 
process with current, (c) a pure dephasing process without 
current flow and (d) a current-carrying dephasing process. 

Figure [5] shows a choice of processes entering the Red- 
field tensor. All processes contribute to dephasing, be- 
cause the phase of an electron, which is injected from 
the reservoirs, is always random. Figures (a) and (b) 
illustrate relaxation processes. Only (b) contributes to 
the current, i.e., in general the relaxation rate must not 
be confused with the cotunneling current. In (c) and 
(d) two pure dephasing processes are presented, only (d) 
contributes to the current flow. In general, processes 
without current can emerge, if the cotunneling processes 
take place between a single lead and the two-state sys- 
tem (TSS). The observable current is then given by the 
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difference of current-carrying processes in forward and 
backward direction. We have evaluated the rates enter- 
ing, Eqs. Qand 111 using the s. Due to the high 

k,k,c,d 

number of terms, details are not shown and will be given 
elsewhere-^. 

We now turn to the discussion of our results, starting 
with the relaxation time r r . We observe in Fig. [3] that 



1x10 



8x10 - 



w 6x1 6 
E 



% 4x10" - 



2x10 



4X10 
2x10 s 




E a 

E a 



=8.6 neV, y=8.6 neV 
=8.6 neV, 7=17.2 (±eV 
=8.6 (ieV, 7=25.8 neV 
=0 (4.eV, y=17.2 |ieV 
=-8.6 neV, 7=17.2 neV 




-1.5 



-1.0 



-0.5 0.0 

Voltage V/28 



0.5 



1.0 



1.5 



FIG. 3: (Color online) Relaxation time r r in units of Tg = 
the period of coherent oscillations between the two 
molecular states. Different values of e as and 7 are taken, when 
the bias voltage V/2S is varied (with /^ av = + /al)/2 = 



75.832 ^eV and fc B T/> av = 1.136 
0.159. 



l(T d ); inset: fc s T/^ av = 



for an asymmetric TSS, i.e., for e as 7^ 0, there is a pro- 
nounced peak of the relaxation time aiV = — sgn(e as )25, 
where V — hr — fJ-L, he. the sign has to be chosen with 
opposite polarity to the asymmetry energy. This means 
in particular that the relaxation is minimal far away from 
equilibrium. This is the central result of our paper. It is 
most clearly visible for T -C 26, but obviously still domi- 
nates the calculated result for temperatures T w 26, as it 
can be seen in the insets of Figures |21 and In order to 
remain in the cotunneling regime, the voltages are still 
quite small as compared to the excitation energy to the 
next charge states e V2 and e Vo , but on the order of the 
molecular level splitting, i.e. |V| « 26 <C e V2 ,e Vo . For 
quantum computation, achieving a maximum relaxation 
time is e.g. appreciable during read-ou&. 

Although surprising, it can be understood from the 
analysis of the different rates, that V = does not neces- 
sarily imply the lowest relaxation rate. At V = there is 
no net current, i.e. no net exchange of particles in the en- 
semble average, however, this is achieved by the cancella- 
tion of finite currents of equal size in forward and back- 
ward direction. These currents are rather small 5 such 
that current heating is reduced to a minimum. To r r , 
Eq. , such current-carrying processes contribute with 
equal sign - the system relaxes no matter to which reser- 
voir. On top of this, one also has to take into account 



the aforementioned current-less relaxation channels. 

The appearance of the peaks as preferred stable points 
in Fig. |3 can be understood based on the analysis of 
the current-carrying processes, [e.g. Figs. [5] (b) and 
(d)] as schematically shown in Fig. |SJ At low voltages, 
\V\ < 26, the system relaxes into a thermal state close to 
the ground state. Relaxation takes place by spontaneous 
emission of energy into the environment and creation of 
an electron-hole pair in the leads. This pair can recom- 
bine through the electrical circuit which fixes the elec- 
trochemical potentials. This leads to electrical current. 
As the voltage is increased away from V = 0, emission 
processes which lead to a current against the polarity of 
the source are suppressed, the others are enhanced, see 
Figs. El (I) and (II). Depending on the asymmetry of the 
double dot, i.e. on the weight of the excited state on the 
left and the right dot, this leads to an enhancement or a 
suppression of the rate. At |V| > 26, the emission pro- 
cesses against the source are completely blocked: the dot 
relaxation does not provide enough energy to overcome 
the electromotive force. The rate vanishes linearily as 
a function of voltage reflecting the size of the available 
phase space for cotunneling, sec Fig. |5] 

At higher voltages, \V\ > 26, inelastic cotunneling 1 ' 
sets in, see Figs. |S] (III) and (IV): The source provides 
enough energy to even excite the double dot, creating a 
nonequilibrium steady-state population of the molecular 
levels. Hence, inelastic cotunneling provides a way for 
the dot to absorb energy from the environment even at 
low temperature. This process can be experimentally 
identified by a sharp increase of the currents-it 

Hence, at V — ±26, three of the four processes de- 
picted in Fig. vanish at low temperatures, whereas at 
V = only two vanish. The linear voltage-dependence of 
the rates leads to the rather sharp cusps seen in Fig. [3] 
This behavior is smeared out at higher temperatures by 
thermal fluctuations. The peak height is set by the re- 
maining processes: Energy emission with the source and 
currentless relaxation, Fig. [21 (a). As explained above, 
the relative weight of the former strongly depends on the 
weight of the excited molecular state on the individual 
dots and thus is responsible for the strong asymmetry of 
the peaks in Fig. |3 for different asymmetry energies. 

Finally, we analyze the properties of the dephasing 
time 70 as a function of the bias voltage. The total de- 
phasing rate contains relaxing as well as flipless ("elas- 
tic") processes. We hence observe in Fig. 01 a peak struc- 
ture at V = — sgn(e as )2<5 as in the relaxation time, Fig. O 
and a similar peak at V = 0. The latter can be under- 
stood from the suppression of flipless processes (energy 
exchange 0)in an analogous way to the relaxation peak in 
Fig. |21 (energy exchange 26) . At low asymmetry energy 
e as < 7, the dephasing time at V — is longest. At high 
asymmetry e as > 7 and at the nonequilibrium working 
point V = — sgn(e as )2<5, is even longer. In general, this 
indicates the existence of two preferable working points 
for quantum computation: One in equilibrium, the other 
again far from equilibrium. As also already seen in the 
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FIG. 4: (Color online) Dephasing time in the same units as 
in the previous figure for different values of e as and 7, when 
the bias voltage V/28 is varied (with /x av = + /Ui)/2 = 



75.832 ^eV and k B T/^ v = 1.136 
0.159. 
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inset of Fig. |3 the voltage dependence at higher temper- 
ature is here smeared out and the peaks merge. 
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FIG. 5: Qualitative voltage dependence of the rates of emis- 
sion [(I), (II)] and absorption [(III), (IV)] processes, see text; 
these rates do not correspond to the processes in Fig. [5] 



A measurement of the relaxation and dephasing times 
should be feasible either by a time-resolved measurement 
of (a z (t)), e.g., through a single-electron transistor or 
point contact, the saturation broadening methodic or 
resonance schemes such as proposed in Ref. |20] for spins. 

Note that parts of the double-dot literature focus 
on decoherence through phonons or photons (see Refs. 
I^ l2lll22 : .23) . whereas we focus on the cotunneling, which 
becomes relevant when phonons are suppressed by a 



1 .5 cavity . If the spin in a dot is used as qubi 
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cotun- 
neling serves as an indirect contribution to decoherence. 

To conclude, we have studied the decoherence of charge 
states in a double quantum dot due to cotunneling. We 
have shown that decoherence can be controlled through a 
bias voltage V (and thus creating a nonequilibrium situa- 
tion) between the two fermionic baths. In particular, the 
optimum working point for read-out and potentially also 
for operation of the qubit can be in an out-of-equilibrium 
situation at a voltage V = — sgn(e as )2(5. We have given a 
consistent physical interpretation of our findings in terms 
of stability and phase space. This effect of stabilization 
through nonequilebrium should potentially be significant 
for other qubit candidates as well. 
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